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We present the first calculation oi Z + A jet production with heavy 
^ ■ quark fiavor identification at the Tevatron pp collider. The Z + 4 

jet channel is especially interesting as a normalizer for the W + 4 
jet background to top quark signals, as a background to a possible 
t ^ cZ fiavor-changing neutral-current (FCNC) decay signal, and as a 
background to missing-p^ signals from gluino pairs. We also calculate 
the contributions to Vr + 4 jet production from all the different heavy- 
flavor final states. The MADGRAPH program is used to generate 
all leading order subprocess helicity amplitudes. We present Monte 
Carlo results with separation and acceptance criteria suitable for the 
Tevatron experimental analyses. The dependence of the cross sections 
on experimental cuts and the theoretical ambiguities due to the scale 
dependence are discussed. The predicted (1^-1-4 jet)/(Z-|- 4 jet) ratio 
is insensitive to most of these choices. 



There are many potential new physics processes at hadron coUiders, that 
would lead to final states with a weak boson plus multi-jets, where the weak 
boson is identified by its leptonic decay; these signals sometimes also contain a 
second weak boson, whose hadronic decay is less easily identified. Since a weak 
boson can also be produced along with gluon and quark jets, a knowledge of 
these QCD backgrounds is essential to the identification of new physics signals. 
Considerable effort has been devoted in recent years to the calculation of QCD 
W + n iet {n = 1, 2, 3, 4) and Z + n jet (n = 1, 2, 3) cross sections; for the cases 
of high jet multiplicities n, that would be given by many interesting new physics 
signals, these calculations can currently be made at tree level only |jl|, 0, ^ ||. 
We present first results for Z production with four QCD jets, evaluated for the 
Tevatron pp collider at ^/s = 1.8 TeV, including a separation of contributions 
from different heavy quark fiavors. We also calculate W + 4 jets with heavy 
quark fiavor identification; this goes beyond previous PF + 4 jet results that fiag 
only 6-fiavor |jl|. 

Major areas of physics interest in a QCD Z+4 jet calculation are the following. 

(a) The most immediate interest is related to the top-quark search at the Tevatron 
1^, 1^, where in the single-lepton signal with a 6-tag the QCD VT -|- 4 jet channel 
gives the major background, and a comparison of the W/Z ratio could provide 
a calibration; this ratio should be insensitive to theoretical uncertainties in the 
individual cross sections. Furthermore, experimental acceptance and detector 
effects are also expected to cancel in the ratio. By calculating separate cross 
sections for different final-state quark fiavors, we are able to apply our results to 
the case where a heavy quark is tagged. 

(b) Possible isosinglet heavy quarks x would have both charged- current and 
neutral- current decay modes, with branching fraction ratios ||^, ||, p, |10|, |lT| 



B{x -^ qW) : B{x -^ q' Z) ~ 2 : 1. (1) 

(c) A related question is the possible existence of a prominent FCNC decay mode 
of the top quark [^, [1^, |I4|, |T5|, t ^> cZ along with the standard t — > hW decay. 
In this scenario, tt production would lead to a tt — > {cZ){hW) — > Z-l-4 jet signal, 
that must be distinguished from QCD background. 

(d) The production of supersymmetric particles gives rise to missing-pr plus 
multijet signals at hadron colliders. In particular, production of gluino pairs 
gg with decays g — * Xil'i to the lightest neutralino x? are expected to be a source 
of missing-pT plus 4 jets. Here Z-l-4 jet production with invisible Z ^* vv decays 
is the dominant standard physics background. In the case of 6-tagged events, 
there are regions of parameter space where g ^ tt or g ^ bb decays are dominant 



g. 



We now turn to the method used in our Z -|- 4 jet calculation. An impedi- 
ment in calculating subprocesses with many final partons is the large number of 
Feynman diagrams to be enumerated and expressed as amplitudes. For example. 



the gg -^ Zqqgg subprocess involves 516 diagrams. This phase of the calculation 
can be accomplished for any given subprocess by the MADGRAPH program []17| , 



which automatically generates all Feynman graphs and their helicity amplitudes, 
employing the HELAS approach |jl8|. However, MADGRAPH does not enumer- 



ate the contributing subprocesses, which must be entered individually, nor does 
it carry through the cross section calculation, folding in initial parton distribu- 
tions and final phase space integration. We have added a phase-space generator 
and folded in the MRS set D'_ parton distributions [0, evaluated at a scale 
Q^ = (pt) + M"^, where (pt) is the average transverse momentum of the par- 
tons. The renormalization scale in as is set equal to Q^ and the A value is chosen 
accordingly to the value in the parton distribution functions with five flavors. A 
similar procedure is followed in our VT -|- 4 jet calculations. 

For semi-realistic simulations, we make parton-level calculations of pp — > 
W{Z) + 4 jets at ^/s = 1.8 TeV. We identify final partons with jets when 

PtU) > 20 GeV, \v{j)\ < 2, AR{jj) > 0.4, (2) 

where [Ai?(jj)]^ = [Ari{jj)]'^ + [A0(jj)]^ defines the angular separation between 
two jets. A correction must be made in comparing the parton transverse momen- 
tum Pt with the observed (uncorrected) jet transverse energy Et; according to 
CDF simulations p, typically 5 GeV or more must be added to the latter. A full 
simulation including fragmentation and detector characteristics must be made for 
detailed comparisons with experiment. 

For the case that Z is detected by Z — > ee and W is detected by VT — >• eu, we 
take the electron and missing transverse momentum -^j. acceptance to be 

PT(e)>20GeV, |r7(e)| < 1, (3) 

^j. > 20 GeV (for W events), (4) 

and require that the electrons are isolated from jets by AR{ej) > 0.4. These ac- 
ceptance criteria approximate but do not exactly duplicate those used in Tevatron 

experimental analyses. 

Unless otherwise stated, in the following Z denotes Z -^ e^e~ and W denotes 
W^ ^ e^z/ ; with these leptonic branching fractions included, the cross sections 
times branching fractions are denoted Ba. Comparison with experiment requires 
the inclusion of instrumental efficiencies also. 

The total cross sections with these acceptance criteria are 

5a(Z + 4jet) =20.5fb, 5cT(Vr + 4jet) = 337 fb. (5) 

These predictions for the absolute cross sections are however sensitive to the 
choice of the scale Q^ in a^. For Q"^ = (px) the cross sections are higher by a 

factor 2.02 (1.93) for Z + 4 jet {W + A jet) production. 

The relative numerical contributions to the total cross section from different 

subprocesses according to the number of quarks involved in the process are (in 

percentages): 
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(6) 

Thus the six-quark contributions are not very important. The percentage con- 
tributions from the different initial state configurations are 



(7) 



The gluon-gluon initiated production is rather insignificant at the Tevatron en- 
ergy. 

Tagging of 6-quarks is an important means of selecting final states such as 
tt and gg^ containing heavy quarks. In the CDF top-quark search, two means 
of tagging are employed, a silicon vertex detector (SVX) and a soft-lepton tag 
(SLT); the former identifies displaced vertices and the latter identifies leptons 
from h — >■ ivX and 6 — ;> c — ;> ivX. However, the situation is complicated by the 
possibility that a c-jet or a light parton jet may be mistagged as a 6-jet. The 
probability of tagging any particular final state therefore depends on the separate 
probabilities e^ that any single jet j = 6, c, qj g satisfies the tagging criteria. In our 
later assessments of tagging, we will assume the values e^ = 0.18 (e.g. 0.11 from 
SVX and 0.07 from SLT), ec = 0.05 and e^/g = 0.01, which are approximately the 
efficiencies in the CDF top-quark search 0. 

For application to tagging studies, we present here the cross sections in fb for 
different final flavor conflgurations. 
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Folding in the 6-tagging efficiencies given above, we obtain the following tagged 



cross sections: 

no. of tags Bo-(Z + 4jets) Ba{W + Aiets) 

> 20.5 337. 

> 1 1.23 18.1 (9) 

> 2 0.057 0.68 

> 3 0.001 0.01 

To compare these numbers with experimental event rates, one has to multiply 
these cross sections by efficiency factors for the electrons and muons and also 
take into account effects of detector simulations. However, these effects (together 
with theoretical uncertainties) are expected to cancel approximately in the ratio 
of {W + 4jet)/(Z + 4jet) cross sections. 

The predicted W/Z ratio in 4-jet events with at least one 6-tag is about 14.7. 
In contrast to the individual cross sections, this ratio is fairly insensitive to the 
choice of the scale Q^ in a^. With Q^ = (pt) the ratio is 14.1. This number is 
also fairly insensitive to the jet threshold pt cut. Even if we relax the pt and r) 
requirements on the fourth jet (as CDF does to increase statistics in the top-quark 
sample), this W/Z ratio remains about 14. This ratio does however depend on 
the lepton rapidity cut; for \ri{£)\ < 2.5 we obtain a W/Z ratio of about 10. 

With 19.2 pb~^ luminosity, CDF finds two 6-tagged Z + 4 jet events and 
seven 6-tagged ly + 4 jet events, with relaxed Et and r] requirements on the 
fourth jet. Although the statistics are small, this observed W/Z ratio in 4-jet 
events appears to be anomalously low in comparison with the QCD prediction. If 
future statistics conffim that the 6-tagged WAj/ZAj ratio is indeed significantly 
lower than the pure QCD ratio, then there must be new physics in the Z -f 4-jet 
channel. Furthermore, if the tagged W + 4-jet events are indeed dominated by ti 
production, as suggested by the CDF analysis ^ and by our results above, then 
the tagged Z + 4-jet eventsQ are dominated by new physics beyond the standard 
model. 

Interesting possibilities for such new physics include (i) a singlet charge -1/3 
quark Xb, that mixes with the 6-quark and therefore has a prominent Xb — >■ bZ 
decay mode P, or (ii) FCNC decays of the top quark t — > cZ that would follow 
from mixing of t with a charge 2/3 singlet quark [Q. In case (ii), the 6-tag would 
have to be faked by the c-jet. 

We next consider the QCD {Z -^ vv) + 4 jet background to the missing-p^ 
signals of supersymmetry. We here consider missing-pr requirements of 

i)^ > 50 or 100 GeV , (10) 

along with the same jet cuts as before. The integrated cross sections are 

^ The contribution to the Z + 4 jet final state from pp ^ W + Z + 2 jets ^^ with hadronic 
W decay is much smaller than the number from the Z + 4 jet channel. 



Ba{i)^ > 50; 4 jets) = 283 fb Ba{:^^ > 100; 4 jets) = 97 fb (11) 

The contribution for different final flavor conflgurations are: 

b c q/g 5o-(|^y > 50 + 4jets) 5(T(|^y > 100 + 4jets) 
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Including the tagging efficiences assumed above, the tagged cross sections are as 
follows. 

no. of tags Ba{^r^ > 50 + 4 jets) Ba^^r^ > 100 + 4 jets) 

> 283. 97. 

> 1 16.1 5.56 (13) 

> 2 0.70 0.24 

> 3 0.01 0.004 

A detailed consideration of the dynamical distributions of Z + 4-jet events will 
be presented elsewhere. 
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